Journal  of  Power  Sources  253  (2014)  123-131 


ELSEVIER 


Contents  lists  available  at  ScienceDirect 

Journal  of  Power  Sources 

journal  homepage:  www.elsevier.com/locate/jpowsour 


Three-dimensional  study  of  compressed  gas  diffusion  layers 
using  synchrotron  X-ray  imaging 

C.  T6tzkea,b’*,  G.  Gaiselmannc,  M.  0senbergb,  J.  Bohnerd,  T.  Arlt  ,  H.  Markotter  , 
A.  Hilger  ,  F.  Wieder3'  ,  A.  Kupsch6,  B.R.  Muller e,  M.P.  Hentschel e,  J.  Banhart3'  , 
V.  Schmidt c,  W.  Lehnertd' ,  I.  Manke 

aTechnische  Universitdt  Berlin,  10623  Berlin,  Germany 

b  Helmholtz-Zentrum  Berlin  fiir  Materialien  und  Energie,  14109  Berlin,  Germany 
c  Universitdt  Ulm,  89069  Ulm,  Germany 
d  Forschungszentrum  Jtilich,  52425  Jiilich,  Germany 

e  BAM  Bundesanstalt  fiir  Materialforschung  und  -priifung,  12200  Berlin,  Germany 
fRWTH  Aachen  University,  52056  Aachen,  Germany 


CrossMark 


HIGHLIGHTS 


•  Microstructure  and  shape  of  GDL  were  studied  using  synchrotron  X-ray  tomography. 

•  Geometric  properties  alter  significantly  upon  compression. 

•  Direct  impact  on  transport  conditions  for  both  gas  and  fluid  flow  in  the  GDL. 

•  Tortuosity  and  connectivity  of  3D  transport  paths  through  the  pore  phase  were  studied. 

•  Pronounced  fiber  extension  into  the  flow-field  channel  volume  upon  compression. 
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We  present  a  synchrotron  X-ray  tomographic  study  on  the  morphology  of  carbon  fiber-based  gas 
diffusion  layer  (GDL)  material  under  compression.  A  dedicated  compression  device  is  used  to  provide 
well-defined  compression  conditions.  A  flat  compression  punch  is  employed  to  study  the  fiber  geometry 
at  different  degrees  of  compression.  Transport  relevant  geometrical  parameters  such  as  porosity,  pore 
size  and  tortuosity  distributions  are  calculated.  The  geometric  properties  notably  change  upon 
compression  which  has  direct  impact  on  transport  conditions  for  gas  and  fluid  flow.  The  availability  of 
broad  3D  paths,  which  are  most  important  for  the  transport  of  liquid  water  from  the  catalyst  layer 
through  the  GDL,  is  markedly  reduced  after  compression.  In  a  second  experiment,  we  study  the  influence 
of  the  channel-land-pattern  of  the  flow-field  on  shape  and  microstructure  of  the  GDL.  A  flow-field 
compression  punch  is  employed  to  reproduce  the  inhomogeneous  compression  conditions  found  dur¬ 
ing  fuel  cell  assembly.  While  homogenously  compressed  underneath  the  land  the  GDL  is  much  less  and 
inhomogeneously  compressed  under  the  channel.  The  GDL  material  extends  far  into  the  channel  volume 
where  it  can  considerably  influence  gas  and  fluid  flow.  Loose  fiber  endings  penetrate  deeply  into  the 
channel  and  form  obstacles  for  the  discharge  of  liquid  water  droplets. 

©  2013  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

Gas  diffusion  layers  (GDLs)  play  a  key  role  for  the  efficient 
operation  of  fuel  cells.  The  water  management  of  proton-exchange 
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fuel  cells  (PEFCs)  crucially  relies  on  both  water  and  media  trans¬ 
port  inside  the  gas  diffusion  layers.  Common  GDL  types  consist  of 
porous  carbon  fiber-based  materials  that  allow  for  a  concurrent 
transfer  of  gaseous  reactants  from  the  flow-field  channels  to  the 
electrodes  and  an  effective  removal  of  liquid  water  produced  at 
the  catalyst  layer.  Furthermore,  the  GDL  must  keep  the  membrane 
humidified  to  maintain  its  proton  conductivity  [1—3].  Therefore, 
an  important  design  requirement  of  GDL  materials  is  a  careful 
adjustment  of  water  transport  and  storage  capacities.  Numerous 
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theoretical  4-8]  and  experimental  [9,10]  studies  focusing  on  two- 
phase  flow  in  porous  media  contributed  to  the  basic  under¬ 
standing  of  underlying  water  transport  mechanisms.  Recently, 
neutron  [11,12]  and  X-ray  imaging  [13—17]  were  shown  to  be  well 
suited  for  the  investigation  of  water  in  porous  materials  such  as 
GDLs  [18-28].  The  transport  properties  of  GDLs  are  closely  linked 
to  their  microstructure.  Therefore,  a  precise  knowledge  of  the 
three-dimensional  morphology  is  a  basic  prerequisite  to  under¬ 
stand  the  functionality  of  GDLs. 

A  sophisticated  way  of  systematically  designing  GDL  materials  is 
to  detect  microstructures  with  improved  physical  properties  by 
means  of  model-based  computer  experiments.  The  development  of 
appropriate  stochastic  microstructure  models  allows  one  to  create 
virtual  3D  structures,  the  transport  properties  of  which  are  subse¬ 
quently  evaluated  by  means  of  numerical  transport  simulations. 
Systematic  parameter  modification  of  the  microstructure  model 
together  with  numerical  transport  simulation  can  help  to  identify 
most  suitable  transport  geometries  [29,30].  The  combination  of  3D 
analysis,  geometric  modeling  [31,32  and  transport  simulation 
[33,34]  is  considered  a  promising  approach  to  establish  a  sound 
basis  for  systematic  virtual  design  of  fiber-based  porous 
microstructures. 

We  performed  high-resolution  synchrotron  X-ray  measure¬ 
ments  on  various  GDL  materials  at  different  degrees  of  compres¬ 
sion.  This  broad  experimental  data  set  serves  as  input  for  the 
development  and  validation  of  stochastic  microstructure  models.  In 
this  paper,  we  focus  on  the  nonwoven  fiber-based  material  H2315 
produced  by  Freudenberg  FFCCT. 

2.  Methods  and  materials 

2.1  Imaging  conditions 

The  experiments  were  performed  at  the  Helmholtz-Zentrum 
Berlin  using  the  synchrotron  tomography  station  of  the  BAMline 
which  is  located  at  the  synchrotron  source  BESSY  II  (Berlin/Ger¬ 
many)  [35].  The  principal  set-up  of  the  tomographic  measurement 
is  displayed  in  Fig.  1. 

A  W— Si  multilayer  monochromator  was  employed  to  turn  the 
white  beam  into  monochromatic  X-ray  light  with  an  energy  reso¬ 
lution  of  A E/E  =  10  2.  The  beam  energy  was  adjusted  to  15  keV  in 
order  to  achieve  optimal  contrast  for  fibers.  A  PCO  camera 
(4008  x  2672  pixel2)  in  combination  with  a  lens  system  and  a 


rotation  table 


Fig.  1.  Scheme  of  the  principal  experimental  set-up  of  the  tomographic  instrument. 


20  pm  thick  CWO  scintillator  screen  was  use  to  capture  a  field  of 
view  of  3.6  x  2.3  mm2  which  corresponds  to  an  image  pixel  size  of 
0.876  pm  and  a  respective  physical  spatial  resolution  of  about  2  pm 
[36].  Circular  GDL  samples  of  3  mm  diameter  were  aligned  by  the 
sample  holder  and  mounted  on  a  translation/rotation  stage.  For 
each  tomogram,  a  sample  was  rotated  stepwise  over  an  angular 
range  of  180°.  A  radiographic  set  of  1500  projections  and  500  flat- 
field  images  was  acquired  and  subsequently  reconstructed  to  a  3D 
image  volume.  The  exposure  time  was  2.5  s  plus  1.7  s  read-out  time 
for  a  single  radiograph  adding  up  to  an  acquisition  time  of  140  min 
for  an  entire  tomogram. 

2.2.  Sample  compression  device 

Fligh-resolution  synchrotron  studies  on  compressed  micro¬ 
structures  of  fuel  cell  components  such  as  GDL  or  membrane 
electrode  assemblies  call  for  a  dedicated  compression  device  that 
ensures  well-defined  compression  conditions  during  tomographic 
measurements.  The  design  of  this  device  must  allow  for  both 
planar  and  horizontal  alignment  of  a  sample  and  precise  adjust¬ 
ment  of  its  thickness  with  an  accuracy  of  a  few  micrometers. 
When  aligned  in  planar,  horizontal  manner,  the  sample  orienta¬ 
tion  coincides  with  the  reconstructed  tomographic  slices  which 
has  practical  advantages  for  the  geometric  analysis.  In  order  to 
create  different  compression  conditions  the  device  is  equipped 
with  differently  shaped  compression  punches:  A  flat  punch  re¬ 
alizes  homogenous  compression  while  a  punch  with  an  integrated 
channel  profile  simulates  the  inhomogeneous  compression  exer¬ 
ted  by  the  flow-field  during  fuel  cell  assembly  which  gives  rise  to 
an  embossed  channel-rip-pattern.  Moreover,  some  instrumental 
requirements  have  to  be  taken  into  account:  sufficient  beam 
transmittance  has  to  be  ensured,  dimensional  restrictions  with 
respect  to  the  small  field  of  view  apply  (usually  some  millimeters 
at  a  spatial  resolution  of  about  1  pm),  and  mechanical  stability  is 
required  in  order  to  keep  the  sample  at  rest  during  tomographic 
scans. 

In  order  to  satisfy  these  requirements  a  dedicated  compression 
device  has  been  designed  and  constructed  at  the  For- 
schungszentrum  Jiilich.  It  facilitates  tomographic  ex  situ  measure¬ 
ments  of  GDL  components,  such  as  GDL  or  membrane  electrolyte 
assemblies  at  the  synchrotron.  Due  to  its  chemical  insensitivity,  e.g. 
to  acids,  the  device  can  be  also  used  to  study  components  of  high- 
temperature  PEM  fuel  cells  which  contain  phosphoric  acid  [37—39]. 

Fig.  2  illustrates  the  principal  design  of  the  compression  device. 
The  base  of  the  device  is  depicted  in  Fig.  2B  (left).  It  features  a 
platform  on  which  the  circular  sample  is  placed  on.  The  main 
compression  unit  including  the  adjustable  compression  tool  is 
fastened  to  the  top  of  the  base  element.  This  design  facilitates  quick 
and  easy  exchange  of  samples  which  helps  to  save  valuable 
beamtime  at  the  electron  storage  ring.  An  adjusting  screw  with  an 
ultra-fine-pitch  thread  allows  for  a  vertical  adjustment  of  the 
compression  punch  with  a  precision  of  about  ±5  pm.  The 
compression  punch  does  not  twist  during  vertical  adjustment,  thus 
providing  a  torsion-free  sample  compression.  Relative  vertical 
displacements  can  be  read  from  a  scale  visually,  while  the  absolute 
position  of  the  punch  is  verified  by  analyzing  the  radiographs. 
Within  the  field  of  view  the  device  is  made  of  the  durable  high- 
performance  polyimide-based  polymer  Vespel®  with  a  reduced 
wall  thickness  of  1  mm  which  ensures  high  beam  transmission. 

The  compression  unit  is  equipped  with  exchangeable 
compression  tools:  a  flat  punch  generates  homogenous  compres¬ 
sion  and  a  flow-field  punch  is  used  to  simulate  inhomogeneous 
compression  conditions  as  found  in  assembled  fuel  cells.  The  latter 
punch  features  a  0.8  mm  wide  and  1  mm  deep  channel  profile  (see 
Fig.  2B  right). 
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Fig.  2.  Compression  device  used  for  high-resolution  tomography  of  circular  GDL  samples.  A:  perspective  view  showing  the  interior  construction  of  the  device.  B:  Sample  holder 
basis  (left)  and  compression  cylinder  with  integrated  flow-field  channel  geometry  (right).  C:  Detail  of  the  sample  bearing  section. 


2.3.  Materials 

The  present  tomographic  study  covers  a  range  of  different  GDL 
materials  measured  at  four  stages  of  increasing  compression.  In  this 
paper,  we  focus  on  the  GDL  type  H2315  which  is  a  nonwoven, 
carbon  fiber-based  felt  material  produced  by  the  company  Freu- 
denberg  FFCCT.  The  diameter  of  carbon  fibers  used  in  this  material 
is  about  10  pm.  This  basis  material  is  used  as  a  standard  GDL  in 


many  fuel  cells  where  its  physical  properties  are  adapted  to  the 
requirements  of  the  respective  application,  e.g.  by  adding  water 
proofing  compounds  and/or  a  microporous  layer  (MPL). 

3.  Results  and  discussion 

During  fuel  cell  assembly  GDLs  are  typically  compressed  by 
about  20-30  vol.%.  The  actual  degree  of  compression  depends  on 


Fig.  3.  3D  representation  of  a  nonwoven  GDL  sample  (type  Freudenberg  H2315)  at  different  degrees  of  compression  0%,  9%,  21%,  31%  (from  top  to  bottom),  sample  is  rendered  in  an 
area  of  1  x  1  mm2. 
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the  clamping  pressure  and  elasticity  of  the  employed  GDL  material. 
In  this  section  the  microstructure  of  GDL  samples  compressed  by 
both  a  flat  and  by  a  flow-field  punch  is  analyzed  at  four  stages  of 
increasing  compression  (about  0  vol.%,  10  vol.%,  20  vol.%,  and 
30  vol.%). 

3.1.  Homogenous  compression 

In  Fig.  3  the  morphology  of  the  homogenously  compressed  GDL 
sample  is  rendered  in  an  area  of  1000  x  1000  pm2.  Using  the 
compression  device,  the  initial  thickness  of  the  sample  Zo  =  223  pm 
is  successively  reduced  to  Zio  =  203  pm,  Z20  =  176  pm,  Z30  =  154  pm, 
which  corresponds  to  an  actual  degree  of  compression  of  9  vol.%, 
21  vol.%,  31  vol.%,  respectively.  We  refer  to  the  initial  thickness  Zo  as 
“uncompressed  state”  even  though  it  actually  includes  a  pre¬ 
compression  of  about  0.025  MPa,  which  is  necessary  to  create  a 
plain  sample  orientation  and  a  smooth  GDL  surface  which  is 
needed  to  be  able  to  determine  the  initial  GDL  thickness.  As 
apparent  in  Fig.  3  the  tomographic  data  is  sufficiently  well  resolved 
to  clearly  follow  the  course  of  single  fibers  and  to  distinguish  be¬ 
tween  fibers  even  if  they  are  located  very  close  to  each  other. 
A  comparison  of  the  3D  images  shows  that  the  compression  of  the 
GDL  material  takes  place  at  the  expense  of  the  pore  volume.  No 
notable  lateral  dilatation  of  sample  was  observed,  suggesting  that 
fibers  are  almost  exclusively  translated  in  the  z-direction  of 
compression. 

In  order  to  analyze  the  influence  of  compression  on  micro- 
structural  characteristics  we  computed  the  porosity,  pore  size  dis¬ 
tribution  (PSD)  and  geometrical  tortuosity  of  the  material  for 
different  degrees  of  compression.  Fig.  4A  shows  that  the  initial 
porosity  of  about  78%  of  the  uncompressed  material  decreases  to 
67%  at  a  compression  degree  of  31%,  which  expresses  the  fact  that 
only  the  pore  volume  is  reduced  whereas  the  fiber  volume  remains 
constant. 

Another  characteristic  of  porous  materials  is  the  PSD.  In  fiber 
structures,  the  entire  pore  space  is  interconnected  and  cannot  be 
segmented  into  individual  pores  due  to  the  lack  of  a  physical  cri¬ 
terion.  For  this  reason,  we  analyze  the  pore  structure  in  terms  of  the 
continuous  pore  size  distribution  (cPSD),  which  does  not  consider 
pores  as  discrete  objects  but  treats  the  entire  pore  structure  as 
continuum.  The  cPSD  describes  the  pore  volume  P(r)  that  can  be 
covered  by  spheres  with  a  radius  r,  such  that  these  spheres  do  not 
overlap  with  the  solid  (fiber)  phase.  P(r)  is  calculated  as  a 
normalized  volume  with  respect  to  the  entire  volume  and, 


pore  radius  /  |jm 

Fig.  5.  Difference  in  cPSD  between  the  uncompressed  sample  and  the  one  compressed 
by  31%. 


therefore,  lim  P(r)  =  porosity.  For  a  more  detailed  mathematical 
description  of  the  cPSD  we  refer  to  Ref.  [40].  Fig.  4B  demonstrates 
that  the  compression  of  the  GDL  material  results  in  a  shift  of  the 
entire  cPSD  towards  smaller  pore  radii.  Flowever,  even  at  the 
highest  degree  of  compression  applied  more  than  25%  of  the  total 
volume  consists  of  pores  with  a  radius  of  at  least  10  pm.  Pores  larger 
than  30  pm  can  hardly  be  found  for  any  of  the  compression  states. 

In  Fig.  5,  the  difference  between  the  cPSD  for  the  uncompressed 
sample  and  the  one  compressed  by  31  vol.%  is  shown.  Starting  from 
large  pore  radii  the  difference  in  cPSD  increases.  The  maximum 
difference  is  found  at  a  pore  radius  of  about  10  pm,  after  which  the 
difference  decreases  towards  smaller  radii.  This  proves  that  the 
compression  of  the  pore  phase  occurs  mostly  at  the  expense  of 
pores  larger  than  R  =  10  pm.  Naturally,  smaller  pores  are  also 
affected  by  compression.  Flowever,  this  is  outweighed  by  the  con¬ 
version  of  large  pores  into  smaller  ones.  When  considering  pores 
ranging  from  5  pm  to  10  pm  radius  only,  we  find  that  such  pores 
make  up  31%  of  the  total  volume  in  the  compressed  and  only  25%  in 
the  uncompressed  state. 

Another  important  structural  characteristics  with  respect  to  the 
transport  of  water  and  gas  is  the  tortuosity  of  the  pathways  through 
the  pore  phase.  We  investigated  the  ‘geometric  tortuosity’  of  the 
pore  space  in  through-plane  direction,  which  is  defined  as  the 


Fig.  4.  A:  Porosity  plotted  against  compression  of  the  GDL.  B:  Plot  of  cPSD  for  different  degrees  of  compression.  Values  are  calculated  on  an  area  of  700  x  700  pm2. 
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Table  1 

Geometric  tortuosity  calculated  in  the  through-plane  direction. 


Compression  [%] 

Through-plane 

connectivity 

Mean 

tortuosity 

Standard 

deviation 

0 

1 

1.50 

0.099 

09 

1 

1.46 

0.090 

21 

1 

1.43 

0.075 

31 

1 

1.48 

0.092 

Euclidean  length  of  the  shortest  paths  of  a  geometric  3D  graph 
representing  all  possible  paths  through  the  pore  phase  divided  by 
the  material  thickness  (in  z-direction).  Starting  from  a  randomly 
chosen  location  on  top  of  the  porous  material,  its  geometric  tor¬ 
tuosity  can  be  represented  by  a  probability  distribution  instead  of 
looking  at  the  mean  tortuosity  only,  see  e.g.  Ref.  [32].  The  geometric 
3D  graph  of  the  pore  phase  is  computed  using  the  skeletonization 
algorithm  implemented  in  the  software  Avizo  7.  The  results  are 
given  in  Table  1  and  Fig.  6.  The  plots  of  geometric  tortuosity  dis¬ 
tribution  are  quite  similar  at  all  considered  degrees  of  compression. 
The  mean  tortuosity  ranges  between  1.42  and  1.50  and  has  no 
obvious  correlation  to  the  degree  of  compression.  To  better  un¬ 
derstand  this  observation  it  is  useful  to  recall  that  the  geometric 
evaluation  of  paths  does  not  yet  include  a  bottleneck  criterion,  i.e. 
the  influence  of  minimum  cross  section  on  the  transport  capacity  is 
not  taken  into  account.  The  course  of  the  paths  is  determined  by  the 
skeletonization  of  the  pore  phase,  i.e.  the  pathways  follow  the 
central  axes  of  pores.  Three-dimensional  representations  of 
percolation  paths  through  the  uncompressed  and  compressed 
material  (see  Fig.  6D)  reveal  that  the  straightness  of  the  paths  is 
very  similar  in  either  case.  Fig.  6B  illustrates  how  a  typical  path  has 
shrunk  upon  compression.  As  visible  in  Fig.  6C  an  original  path 
might  also  be  blocked  by  fibers  displaced  by  compression,  thus 
forcing  an  alternative  transit  pathway  through  the  pore  phase. 

Taking  into  account  the  negligible  transport  capacity  for  liquid 
water  of  very  narrow  transport  paths  we  extend  the  calculation  of 
geometric  tortuosity  by  introducing  a  bottleneck  criterion.  A  min¬ 
imum  lateral  distance  Rmin  of  the  central  axis  of  the  pore  and  solid 
phase  is  required  all  along  the  pathway  in  order  to  be  accepted  as  a 
path.  In  Fig.  7,  the  calculated  mean  geometric  tortuosity  in  through- 
plane  direction  is  displayed  for  three  Rmin:  1,  5  and  10  pm.  The 
diagram  also  includes  the  through-plane  connectivity  of  the  paths, 
i.e.  the  probability  that  starting  from  a  randomly  chosen  point  at 


the  surface  a  continuous  pathway  through  the  GDL  exists.  In  other 
words,  we  determine  the  shortest  passage  of  a  solid  sphere  of 
radius  Rmin  through  the  GDL  via  the  pore  phase.  The  through-plane 
connectivity  is  the  probability  of  a  possible  passage  starting  at  a 
randomly  chosen  point  of  the  surface.  Fig.  7  A  displays  the  results  for 
the  uncompressed  material.  It  is  obvious  that  the  mean  tortuosity 
increases  with  growing  Rmin.  If  a  minimum  radius  of  10  pm  is 
assumed  to  be  the  smallest  constriction  along  the  3D  path  through 
the  GDL  the  mean  tortuosity  amounts  to  about  2.1,  which  is  notably 
larger  than  the  mean  geometrical  tortuosity  of  around  1.5  calcu¬ 
lated  for  Rmin  =  1  pm.  Furthermore,  the  plot  of  the  through-plane 
connectivity  reveals  that  liquid  water  transport  through  most 
pathways  is  limited  by  bottlenecks  with  a  radius  between  5  and 
10  pm.  While  the  probability  of  a  possible  transfer  pathway  for  a 
solid  sphere  with  radius  Rmin  =  5  pm  is  still  about  98%  this  prob¬ 
ability  drops  to  less  than  6%  for  spheres  with  Rmin  =  10  pm.  Fig.  7B 
displays  the  results  for  the  material  compressed  by  31  vol.%. 
Compression  of  the  GDL  aggravates  transport  limitations  by  bot¬ 
tlenecks:  The  through-plane  connectivity  decreases  from  0.98  to 
0.90  for  Rmin  =  5  pm.  For  a  radius  constriction  of  Rmin  =  10  pm  the 
connectivity  even  drops  from  0.056  to  0.014.  When  comparing  the 
mean  tortuosity  for  Rmin  =  10  pm  in  compressed  and  uncompressed 
states  it  appears  that  especially  more  tortuous  paths  are  affected  by 
the  formation  of  bottlenecks  with  a  radius  <10  pm. 

The  existence  of  a  few  wide  paths  is  of  great  importance  for  the 
global  liquid  water  transport.  Due  to  their  low  transport  resistance 
these  paths  provide  preferred  transfer  routes  for  the  removal  of 
liquid  water  from  the  catalyst  layer  [8,41  ].  In  this  study  we  focus  on 
the  through-plane  tortuosity  of  the  material.  This  property  can  be 
regarded  as  a  material  characteristic.  The  in-plane  transport  of 
water  and  gases  in  a  cell  is  also  important,  especially  under  the 
lands.  Due  to  the  anisotropic  nature  of  the  GDL  material,  the  in¬ 
plane  tortuosity  would  differ. 

3.2.  Inhomogeneous  GDL  compression  by  the  flow  field  pattern 

When  assembled  in  fuel  cells,  GDLs  are  inhomogeneously 
loaded  and  compressed  by  the  channel-land-pattern  of  the  flow- 
field.  In  this  section,  we  adapt  the  sample  compression  to  more 
realistic  conditions  using  a  compression  punch  with  an  integrated 
channel— land  profile.  Fig.  8  shows  perspective  views  on  the  GDL 
sample  which  has  successively  been  compressed.  The  initial 


Fig.  6.  A:  Distribution  of  geometric  tortuosity  for  different  degrees  of  compression;  B-D:  3D  representation  of  the  percolation  pathways  through  the  pore  phase  for  compression 
degrees  of  0%  and  31%.  Compressed  paths  are  plotted  in  red,  uncompressed  in  white.  (For  interpretation  of  the  references  to  color  in  this  figure  legend,  the  reader  is  referred  to  the 
web  version  of  this  article.) 
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Fig.  7.  Mean  geometric  tortuosity  at  different  radius  limitations  Rmin  of  percolation  paths  (black  rectangles)  and  the  respective  through-plane  connectivity  (red  circles)  A:  degree  of 
compression  =  0  vol.%;  B:  degree  of  compression  =  31  vol.%.  (For  interpretation  of  the  references  to  color  in  this  figure  legend,  the  reader  is  referred  to  the  web  version  of  this 
article.) 


thickness  Zo  =  219  pm  has  been  reduced  to  Z\q  =  197  pm, 
Z2 o  =  180  pm,  Z30  =  155  pm  which  refers  to  actual  compression 
degrees  of  10  vol.%,  18  vol.%,  and  29  vol.%.  The  reconstructed  image 
volume  was  cut  to  provide  a  cross  sectional  insight.  In  this 
perspective,  it  is  clearly  visible  how  the  contour  of  the  compressing 
tool  is  stamped  into  the  GDL  with  increasing  compression,  even¬ 
tually  generating  considerable  differences  in  the  local  material 
density.  While  the  material  underneath  the  land  is  homogenously 
compressed  to  the  thickness  determined  by  the  compression  de¬ 
vice  the  fiber  structure  underneath  the  channel  is  much  less  and 
inhomogeneously  compressed.  The  material  thickness  along  the 
center  line  of  the  channel  remains  almost  constant.  This  implies 
that  upon  assembly  in  fuel  cells  the  fiber  material  significantly 
extends  into  the  channel  volume  and  occupies  space  at  the  expense 
of  the  media  flow  capacity  of  the  channel.  Fig.  8  also  exhibits  some 
individual  fiber  endings  that  are  bent  up  by  the  channel  edges. 


Fig.  8.  Nonwoven  GDL  sample  compressed  by  a  channel-land  profile.  Compression 
underneath  land  area  is  (from  top  to  bottom)  0  vol.%,  10  vol.%,  18  vol.%,  and  29  vol.%. 


These  fiber  endings  interfere  with  the  gas  flow  inside  the  channel 
and  could  even  alter  the  flow  regime,  e.g.  by  producing  vortexes  or 
by  hampering  the  passage  of  water  droplets.  As  those  fibers  can 
hold  up  moving  droplets  they  are  potential  accumulation  points  for 
liquid  water  in  the  channel. 

We  investigated  the  influence  of  the  flow-field  geometry  on  the 
local  GDL  morphology  by  calculating  geometrical  parameters 
relevant  for  transport  for  regions  located  underneath  both  channels 
and  land.  In  the  following  we  refer  to  these  regions  as  “land  region” 
and  “channel  region”,  respectively.  Fig.  9  illustrates  the  size  and 
position  of  the  analyzed  GDL  sub-regions.  The  dimensions  of  the 
two  land  regions  are  219  x  1314  x  z  pm3,  respectively,  where 
parameter  z  denotes  the  thickness,  which  is  adjusted  to  the 
increasing  compression:  z  is  192  pm,  171  pm,  158  pm  and  131  pm.  In 
contrast  to  land  areas,  the  local  thickness  of  the  GDL  material  in  the 
channel  region  strongly  depends  on  the  distance  to  the  channel 
edge.  The  maximum  thickness  is  found  at  the  channel  center  from 
where  it  continuously  decreases  towards  the  land  area.  With 
increasing  sample  compression,  the  thickness  at  the  center  remains 
almost  constant  or  even  slightly  increases  (see  below  for  more 
details),  whereas  fiber  material  closer  to  the  channel  edges  is 
increasingly  compressed.  We  analyzed  geometric  properties  of  the 
channel  region  within  a  sub-volume  of  701  x  1314  x  245  pm3 
(x  x  y  x  z)  which  was  thoroughly  filled  with  GDL  material  for  all  the 
compression  stages  under  study. 

The  porosity  of  the  GDL  within  both  land  and  channel  region  at 
different  degrees  of  compression  is  displayed  in  Fig.  10.  The  dif¬ 
ference  in  porosity  between  the  two  areas  in  the  uncompressed 
state  is  due  to  the  precompression  of  the  sample  during  prepara¬ 
tion  as  explained  in  Section  3.1.  The  porosity  underneath  the  land 
decreases  according  to  the  reduction  of  the  pore  volume  by  the 
compression  punch  very  similar  to  the  results  of  the  flat  punch 
shown  in  Fig.  4a.  In  contrast,  underneath  the  flow-field  channel, 
were  the  GDL  material  is  not  in  direct  contact  with  the  compression 
punch,  the  porosity  decreases  only  slightly.  This  is  in  good  accor¬ 
dance  with  the  visual  impression  of  the  morphology  provided  by 
Fig.  8. 

In  order  to  study  the  pore  phase  in  more  detail  the  continuous 
PSD  is  calculated  for  both  the  land  and  channel  regions.  A  plot  of 
the  cPSD  is  presented  in  Fig.  11.  Underneath  the  land  the  graphs  of 
cPSD  are  equivalent  to  those  calculated  for  the  homogenous 
compression  by  the  flat  punch.  Comparison  of  the  results  for 
channel  and  land  regions  shows  that  the  channel  region  contains  a 
notable  amount  of  pores  with  a  radius  >30  pm,  whereas  no  pores  of 
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Fig.  9.  Visualization  of  GDL  sub- volumes  analyzed:  Illustration  of  size  and  position 
within  the  sample  (degree  of  compression  29  vol.%).  The  land  regions  are  highlighted 
in  red,  the  channel  region  in  blue.  (For  interpretation  of  the  references  to  color  in  this 
figure  legend,  the  reader  is  referred  to  the  web  version  of  this  article.) 


this  size  are  present  underneath  the  land  area.  When  focusing  on 
radii  >30  pm  in  Fig.  11 A  we  find  that  the  pore  volume  inside  the 
channel  region  slightly  increases  with  compression.  This  seemingly 
paradox  can  be  explained  by  the  concave  deformation  of  fiber 
material  underneath  the  channel.  The  bending  stress  exerted  at  the 
channel  edges  slightly  bows  the  GDL  material.  Consequently,  pores 
located  in  the  top  layers  of  the  sample  are  stretched  while  those 
located  in  the  bottom  layers  are  slightly  compressed  which  results 
in  an  almost  constant  overall  porosity  of  the  channel  region 
considered. 

Fig.  12  shows  the  influence  of  inhomogeneous  compression  on 
the  local  GDL  thickness  and  the  partial  occupation  of  the  flow- 
field  channel  by  the  fiber  material.  The  plot  of  the  sample 
thickness  measured  at  the  channel  center  and  underneath  the 
land  at  different  degrees  of  compression  reveals  that  the  GDL  is 
much  less  compressed  in  the  channel  region.  While  the  thickness 
beneath  the  land  is  successively  reduced  by  10, 18  and  29%  by  the 
compression  punch,  it  even  increases  slightly  at  the  center  line  of 
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Fig.  10.  Porosity  determined  for  land  (red  triangles)  and  channel  areas  (black  dots)  at 
different  degrees  of  compression.  (For  interpretation  of  the  references  to  color  in  this 
figure  legend,  the  reader  is  referred  to  the  web  version  of  this  article.) 
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the  flow-field  channel.  At  29%  compression  the  GDL  is  twice  as 
thick  along  the  center  line  of  the  channel  as  in  the  land  region. 
Consequently,  it  extends  far  into  the  channel  volume  as  displayed 
in  Fig.  12B.  Already  in  the  uncompressed  state  a  penetration 
depth  of  90  pm  is  measured  for  the  GDL.  With  increasing  degree 
of  compression  the  material  continues  to  extend  into  the  chan¬ 
nel.  At  29%  compression  the  penetration  depth  reaches  167  pm. 
Fig.  12C  also  clearly  demonstrates  that  individual  fiber  endings 
protrude  into  the  flow-field  channel.  These  fibers  are  bent  at  the 
edges  of  the  flow-field  channel  profile  and  tilt  up  with  increasing 
degree  of  compression.  At  a  compression  degree  of  29%  some 
fiber  endings  protrude  295  pm  into  the  channel  volume  where 
they  potentially  influence  the  gas  flow  and  form  obstacles  for 
water  droplets  that  are  transported  through  the  flow  field 
channel. 

4.  Conclusions 

Synchrotron  X-ray  tomography  was  used  to  study  the  micro¬ 
structure  of  carbon  fiber-based  GDL  materials  at  different  degrees 
of  compression.  A  dedicated  compression  device  was  used  to  pro¬ 
vide  realistic  compression  conditions  that  reflect  conditions  in  real 
fuel  cells.  The  application  of  different  compression  punches 
allowed  for  studying  geometric  properties  relevant  for  the  trans¬ 
port  of  media  in  homogenously  compressed  GDL  material  as  well  as 
to  investigate  the  influence  of  the  land-channel-pattern  of  the  flow 
field  on  the  local  microstructure  of  the  GDL.  The  main  findings  of 
the  tomographic  study  are: 

1.  The  geometric  properties  of  the  fiber-based  material  are 
altered  significantly  upon  compression.  This  has  a  direct 
impact  on  the  transport  conditions  for  both  gas  and  fluid 
flow  in  the  GDL. 

2.  Compression  mainly  affects  the  pores  while  the  fibers  are 
almost  incompressible. 

3.  When  homogenously  compressed,  fibers  are  mainly  trans¬ 
lated  in  the  z-direction. 

4.  Depending  on  the  degree  of  compression  the  pore  volume 
mainly  consist  of  pores  with  a  radius  ranging  between  5  pm 
and  15  pm. 

5.  We  found  only  a  minor  influence  of  compression  on  the 
geometric  tortuosity  in  through-plane  direction. 

6.  Taking  into  account  the  transport  limiting  effect  of  bottle¬ 
necks  (radius  constrictions  along  the  3D  transport  paths) 
the  mean  tortuosity  of  the  uncompressed  material  in¬ 
creases  from  about  1.5  to  about  2.1  when  only  3D  paths  with 
a  minimal  radius  of  at  least  10  pm  are  included  in  the 
calculation.  Only  a  few  percent  of  the  possible  transport 
paths  through  the  GDL  fulfill  this  criterion  but  that  very 
paths  might  act  as  preferred  water  transport  channels  from 
the  catalyst  layer  to  the  channel.  This  result  is  in  agreement 
with  recent  findings  from  in-operando  studies  by  Markotter 
et  al.  [8,41]. 

7.  The  through-plane  connectivity  of  transport  paths  that 
satisfy  the  bottleneck  criterion  strongly  depends  on  the 
minimum  pore  radius  as  well  as  on  the  degree  of  compres¬ 
sion.  A  compression  of  31  vol.%  decreases  the  number  of 
available  broad  3D  transport  paths  with  a  minimum  diam¬ 
eter  of  at  least  20  pm  to  25%  of  the  initial  value  (see  Fig.  7C). 
These  broad  paths  are  most  important  for  liquid  water 
transport.  Therefore,  compression  causes  a  significant 
decrease  of  low  resistance  water  transport  paths. 

8.  Compression  of  the  fiber  structure  by  the  flow-field 
strongly  alters  the  overall  GDL  structure.  While  homo¬ 
geneously  compressed  underneath  the  land  it  is  much 
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Fig.  11.  Continuous  pore  size  distribution  computed  for  channel  (A)  and  land  regions  (B). 


less  and  inhomogeneously  compressed  underneath  a 
flow-field  channel.  At  the  center  of  flow  field  channels 
even  a  local  increase  in  porosity  might  occur  due  to  the 
pronounced  extension  of  fiber  layers  into  the  channel 
volume. 

9.  The  fiber  material  extends  far  into  the  flow-field  channel 
volume  and  occupies  channel  volume  that  is  actually  pro¬ 
vided  for  gas  flow  and  removal  of  liquid  water.  For  example,  a 


penetration  depth  of  167  pm  was  measured  for  a  compres¬ 
sion  of  29  vol.%. 

10.  We  found  many  fibers  tilted  up  by  bending  stresses  at 
channel  edges  with  their  endings  penetrating  deeply  into  the 
flow  field  channel.  For  example,  a  penetration  depth  of 
295  pm  was  measured  for  a  compression  of  29  vol.%.  This  can 
have  significant  influence  on  both  the  gas  and  liquid  flow  as 
those  fibers  form  obstacles  for  moving  water  droplets  and 
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Fig.  12.  Evaluation  of  GDL  thickness  underneath  land  and  channel  areas.  A)  Plot  of  GDL  thickness  underneath  land  (circles)  and  channel  (triangles).  B)  Plot  of  the  penetration  depth 
of  GDL  (triangles)  and  outermost  fibers  (circles)  into  the  channel  volume.  C)  Perspective  side-view  illustrating  the  thickness  of  the  GDL  sample  at  a  compression  degree  of  29%. 
a  =  thickness  underneath  land;  b  =  thickness  underneath  channel;  c  =  thickness  including  the  most  distant  fiber  ends;  b  -  a  =  penetration  depth  of  GDL  into  channel; 
c  -  a  =  penetration  depth  of  outermost  fiber  ends. 
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can  potentially  create  agglomerations  of  liquid  water  in  the 
fuel  cell.  This  effect  has  also  to  be  taken  into  account  in  future 
flow  field  simulations. 


5.  Outlook 

In  this  study,  the  superior  brilliance  of  a  synchrotron  X-ray  beam 
and  the  high  spatial  resolution  of  the  associated  imaging  instru¬ 
ment  were  exploited  to  generate  tomographic  data  that  reproduce 
even  smallest  morphological  details  of  the  fiber  material  with  great 
accuracy.  This  information  was  used  to  estimate  some  effects  of  the 
microstructure  on  both  gas  and  liquid  transport.  However,  a  more 
detailed  analysis  of  water  and  gas  transport  properties  of  the  GDL 
material  additionally  requires  numerical  transport  computations, 
e.g.  by  Lattice-Boltzmann  simulations.  The  accuracy  of  such  trans¬ 
port  simulations  significantly  relies  on  the  quality  of  geometric 
structure  models  and  the  implementation  of  realistic  boundary 
conditions.  When  simulating  water  and  gas  transport  through  the 
GDL,  the  geometry  of  the  inhomogeneously  compressed  gas 
diffusion  layer  should  be  carefully  taken  into  account.  This  is 
especially  important  when  it  comes  to  computing  mass  and  energy 
transfer  from  the  GDL  into  the  flow-field  channel.  The  highly 
resolved  tomographic  data  produced  in  this  study  is  considered  a 
valuable  source  of  information  for  defining  appropriate  boundary 
conditions  and  is  ideally  suited  as  experimental  basis  for  the 
development  of  stochastic  microstructure  models  needed  for  the 
design  of  materials  with  improved  functionality  by  means  of 
computer  experiments. 
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